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Abstract 

Discrete models have a long tradition in engi- 
neering, including finite state machines, Boolean 
networks, Petri nets, and agent-based models. 
Of particular importance is the question of how 
the model structure constrains its dynamics. 
This paper discusses an algebraic framework to 
study such questions. The systems discussed 
here are given by mappings on an affine space 
over a finite field, whose coordinate functions 
are polynomials. They form a general class of 
models which can represent many discrete model 
types. Assigning to such a system its depen- 
dency graph, that is, the directed graph that 
indicates the variable dependencies, provides a 
mapping from systems to graphs. A basic prop- 
erty of this mapping is derived and used to prove 
that dynamical systems with an acyclic depen- 
dency graph can only have a unique fixed point 
in their phase space and no periodic orbits. This 
result is then applied to a published model of in 
vitro virus competition. 

'This work was supported by the National Science 
Foundation under Grant Nr. CMMI-0908201. 



1 Introduction 

In several areas of engineering, discrete math- 
ematical models, such as finite state machines, 
Boolean networks, Petri nets, or agent-based 
models, play an important role in modeling pro- 
cesses that can be viewed as evolving in discrete 
time, in which state variables have only finitely 
many possible states. Decision processes, electri- 
cal switching networks, or intra-cellular molecu- 
lar networks represent examples. There is a long 
tradition in the engineering literature of work 
related to an understanding of how the struc- 
ture of such models constrains their dynamics, 
see, e.g., [3l|2]. This problem is important both 
for the design of networks as well as for their 
analysis. In both theoretical and applied studies 
of dynamical systems, the problem of predicting 
dynamic features of the system from structural 
properties is very important, see, e.g., [312]. An 
instantiation of this problem is the question in 
systems biology whether the connectivity struc- 
ture of molecular networks, such as metabolic 
or gene regulatory networks, has special features 
that correlate with the type of dynamics these 
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networks exhibit. (This type of application mo- 
tivated the present study.) Some progress has 
been made on this question. For instance, in 
[1] it was shown that certain small network mo- 
tifs appear much more often in such networks 
than would be expected in a random graph. The 
work in this paper was motivated by the desire 
to provide a theoretical framework within which 
the relationship between structure and dynam- 
ics of certain families of dynamical systems can 
be studied. We briefly describe this framework 
and show one example of a result and its impli- 
cations. 

The type of dynamical system studied here 
is discrete in time as well as in variable states. 
That is, we consider a collection x\, . . . ,x n of 
variables, each of which can take on values in a 
finite set X. We consider time-discrete dynami- 
cal systems 

/ : X n — > X n . 

A well-known instantiation of such systems are 
Boolean networks, which have many applica- 
tions in molecular biology as well as engineer- 
ing. In this case, X = {0,1}. As described, 
such systems are set functions without additional 
mathematical structure. It is therefore advanta- 
geous to impose additional mathematical struc- 
ture on X, namely that of a finite field. This 
is of course utilized in the case of Boolean net- 
works. Evaluation of Boolean functions is equiv- 
alent to carrying out arithmetic in the Boolean 
field with two elements. The translation between 
Boolean functions and polynomials is straight- 
forward, with the Boolean operator AND corre- 
sponding to multiplication, addition correspond- 
ing to XOR, and negation corresponding to ad- 
dition of 1. Once we have a finite field structure 
on X, which we shall now denote by K, it is 
well-known that any function K n — > K can be 



represented as a polynomial function [HJ p. 369]. 
Thus, we can focus on systems 

/ = (/i,. ..,/ n ):K»— MK» 

for which each fi is a polynomial function in the 
variables xi, . . . ,x n . We will call such a system 
a polynomial dynamical system (PDS) over K. 
A powerful consequence is that one can use al- 
gorithms and theoretical results from computer 
algebra and computational algebraic geometry, 
such as the theory of Grobner bases. 

It has been shown [TUJ [3] that discrete models 
in several different frameworks can be translated 
into the PDS framework, namely fc-bounded 
Petri nets, so-called logical models used in molec- 
ular biology, and many agent-based models, 
which are becoming very prominent in biology. 
Thus, the framework of PDS provides provides 
a common theoretical formulation. Using al- 
gorithms from polynomial algebra, the software 
package ADAM (Analysis of Dynamic Algebraic 
Models) [6] is very efficient in analyzing various 
features of PDS, including their dynamics. Thus, 
any theoretical results about PDS apply to all 
these model types. 

In this paper, we are concerned with the family 
V of all PDS of a given dimension n over a finite 
field "K. The structural information for a PDS 
includes a directed graph that indicates the de- 
pendency relationships between variables. In the 
context of a molecular network model, this graph 
would represent the wiring diagram of the net- 
work. Graphs can be represented by their adja- 
cency matrices. The set of n x n matrices carries 
an algebraic structure using max and min oper- 
ations, which has been studied previously over 
the real field, see, e.g., Here we prove a re- 
sult about the relationship between properties of 
dynamical systems in V under composition and 
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properties of the corresponding adjacency matri- 
ces in the max-min algebra. As a consequence 
of this result, we derive a dynamic property of 
PDS with acyclic dependency graphs. Finally, 
we give an application of this consequence to a 
published model of in vitro virus competition. 

2 The algebra of dynamical sys- 
tems 

Let K be any finite field. We consider dynamical 
systems 

/ = (A,. ..,/ n ):r->r, 

where fi : K n — > K for i = 1, . . . , n. As observed 
above, any such / is a PDS. Let V denote the 
family of all such systems. It is well-known that 
V has the structure of an associative algebra un- 
der coordinate-wise addition and composition of 
functions. 

To / we can associate a directed graph with 
the n nodes x±, . . . , x n . There is a directed edge 
from Xj to Xi if Xj appears in /j. Let [/] be the 
n x n adjacency matrix of this graph. That is, 
[/] = (dij) is defined as follows: 

f 1 fi depends on Xj 
13 \ otherwise 

Equivalently, cijj = 1 if and only if there are 
p/gel such that 

where x^p) = {x\, . . . , Xj-i,p, x j+i , . . . , x n ) £ 
K n . 

The adjacency matrix has binary entries, and 
we now define two operations on such matri- 
ces. Let B denote the Boolean field with two 
elements, with the natural order < 1. 



Definition 2.1. Given a,b£M, we let: 

1) a©6 = max{a, b} and 

2) a-kb = min{a, b} 

Definition 2.2. Given A = (a^) and B = (bij) 
matrices with entries in B, we define the follow- 
ing operation: 

A * B = C = (dj), where aj = 0fc = i(a ifc * b kj ) 

Remark 2.1. aj = 1 if and only if there is k such 
that aik * b^j = 1 

Definition 2.3. We define A<B if and only if 
dij < bij for all i,j. 

Forming the adjacency matrix of the depen- 
dency graph of a PDS then gives a mapping 

V — ► M. 

And we also have a mapping that associates to 
an element in V its phase space, a directed graph 
on the vertex set |K n |, which encodes the dy- 
namics of the system. We will denote by S the 
set of all directed graphs on |K n | with the prop- 
erty that each vertex has a unique outgoing edge 
(the requirement for being the phase space of a 
deterministic PDS). Hence, we have mappings 

S^V — ► M. 

The result in the next section relates information 
in A4 to information in S. 

3 A property of V — > M 

The main result of this paper is the following the- 
orem. It describes a basic property of the map- 
ping that extracts the adjacency matrix from the 
system. 
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Theorem 3.1. We have that [fog] ^ [/] * [g] 
in M. for all f, g G V . 

Proof. Let A, B, and C be the adjacency matri- 
ces for /, g, and fog, respectively. It suffices to 
prove that, if Cij = 1, then ai k b k j = 1 for some 
k. Let us assume that Cy = 1, then there are 
p 7^ g G K such that 

(fog)i(x^)^(fog)i(xM), 

where = (si, . . . , Xj-i,p,x j+ i, . . . , x n ) G 

K n . Now, since (/ ° g)i = fi o g we have that 

fi(yi, ■■■,yn) fi(yi, ■ ■ ■ ij, 

where y s = g s (x^' p ^) and y s = g s (x^' q ') for s = 
1, . . . , n. Then there is an index k such that ft 
depends on xu and 

Vk ^ Vk, 



or 



Therefore, Ojfc = 1 and bkj = 1 and, conse- 



quently, a ik b kj = 1. 



□ 



Corollary 3.2. [f r ] X [/] r /or a// f € V, and 

for any r > 1. 



Proof. Replace g by / in Theorem 3.1 



□ 



4 PDS with acyclic dependency 
graph 

As an easy application of our main result we now 
consider PDS in V that have an acyclic depen- 
dency graph, that is, no feedback loops in their 
structure. For any / : K n — > K n with acyclic 
dependency graph we can see easily that its ad- 
jacency matrix [/] is a strictly triangular matrix, 



i.e., with zeros in the diagonal (otherwise, there 
would be loops in the graph). Therefore, its 
characteristic polynomial is equal to A r , where 
r is the order of the matrix. So [/] is nilpotent 
and [f] r = 0. 

Corollary 4.1. Any discrete dynamical system 
f : W 1 — > M, n with acyclic dependency graph has 
a unique fixed point. 

Proof. By Corollary [OJ [f r ] ± [f] r = 0, so f r is 
constant: f r = xq. Therefore, xq is the unique 
fixed point of /. 

□ 

Note that the nilpotency index of [/], the 
smallest integer r for which [f] r = and [/] r_1 ^ 
0, gives us an upper bound for the number of 
steps to reach the unique fixed point from any 
other state. 

5 Application 

In this section we will present an application of 



Corollary 4.1 to a published model of an in vitro 
competition between two strains of a murine 
coronavirus studied in Jarrah et. al. The 
model is presented as a hexagonal grid of cells, 
with color coding of cells to indicate their in- 
fection status. Normal cells are represented as 
white. Infected cells are represented as red or 
green, or yellow, in the case of dual infection. 
The infection spreads from the center outwards. 
At each time step one ring of new cells is in- 
fected. The outcome of a cell in the new ring 
depends on the infection status of its two neigh- 
bors in the previous infected ring. The local up- 
date function for each cell is constructed using 
the following rules: 

• If a cell has only one infected neighbor, it 
will get the same type of infection. 
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• If a cell has two infected neighbors, then we 
use the following table to determine the type 
of infection of that cell. 



Rules for the update function 




Green 


White 


Red 


Yellow 


Green 


Green 


Green 


Yellow 


Green 


White 


Green 


White 


Red 


Yellow 


Red 


Yellow 


Red 


Red 


Red 


Yellow 


Green 


Yellow 


Red 


Yellow 



The dynamics of this system is represented in 
Figure [T] In this example we will use 169 cells. 
In order to simulate experimental conditions de- 
scribed in we initialize the infection using 
the 37 center cells. Each cell can have only one 
of four colors at a time. We use the field with 
four elements, F4 = {0, 1, 2, 3}, to represent the 
set of different colors. The color assignment is 
as follows. 



Color assignment 


Color 


Field element 


Green 





Red 


1 


White 


2 


Yellow 


3 



We represent the 169 cells by the variables 
xi, . . . ,xiqq, with xi,...,X37 representing the 
center cells. The variables 2i28> • • • 1^169 repre- 
sents the cells in the outermost ring. 

In |llj the system was studied from the point 
of view of the experimental system. There, a 
collection of cells was infected in the center of 
the dish, and the infection was then observed to 
spread to the rest of the Petri dish in a pattern 
that showed distinct segmentation, matching the 
experimental results. In order to mimic the bio- 
logical system, it was assumed that the cells ini- 
tially infected did not change their infection state 




Figure 1: Dynamics of the constrained model 



subsequently, so that the model is constrained 
and heterogeneous with respect to the rules as- 
signed to all the nodes. This, in effect, changes 
the dynamical system since those cells are now 
assigned constant functions rather than the rules 
described above. The outcome is a steady state 
that shows a characteristic segmentation which 
coincides with experimental observations. See 
Figure [l] for an example. 

If, however, we allow all cells to evolve accord- 
ing to the update rule above, then the outcome 
is significantly different. In that case, the center 
cell plays a very special role. It is represented by 
the only node in the network that does not re- 
ceive any inputs from other cells and is therefore 
constant. (Note that, in order to have an acyclic 
dependency graph, a PDS must have at least one 
node without incoming edges.) The state of that 
node then plays the role of an external parame- 
ter, and a particular choice of state/color for this 
node is part of the system description. In this 
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Figure 2: Adjacency graph for virus competition 



case Corollary |4.1| applies, since the dependency 
graph of the network is acyclic. See Figure [2} 
We conclude that in this case any initialization 
of the system, that is, any assignment of colors 
to each of the nodes in the network, results in a 
unique steady state, namely the state in which 
all nodes have the same color as the center node. 
In effect, what happens is that after initializing 
the 37 nodes in the center they propagate a dis- 
tinct, segmentation-like pattern. At the same 
time, the color in the center cell propagates out- 
ward and overpowers the segmentation pattern, 
resulting in a system state with a homogenous 
color distribution. See Figure [3j 

If we now assume that a larger number of cells 
in the center is infected, that is, is assigned a 
color that does not change subsequently, as in 
the initial experiments, then the color distribu- 
tion along the edge of the infected region propa- 
gates and produces a steady state that show the 
segmentation patterns observed in [TT] . 



The polynomial dynamical system for the 
model in which all cells are allowed to change 
is represented as 



,169 



F. 



169 



/ = (/i,...,/i 69 ) :FJ 

The coordinate functions fa are polynomials in 
F^Xj, Xk], where Xj and x k are the two neigh- 
bors of Xi in the previous infected ring. Let Xj 
represent one of the 169 cells, then X{ is updated 
according fa(xj,x k ), i.e. Xi = fa(xj,x k ), where 
Xj and x k are the two neighbors of Xj in the pre- 
vious infected ring. The following table specifies 
part of the truth table for fa(xj, x&), 



Truth Table 


Color 


Color 


Xj 


Xk 


f (Xj, x k ) 


Color 


Green 


Green 











Green 


Green 


Red 





1 


3 


Yellow 


Green 


White 





2 





Green 


Green 


Yellow 





3 





Green 


Red 


Green 


1 





3 


Yellow 


Red 


Red 


1 


1 


1 


Red 


Red 


White 


1 


2 


1 


Red 


Red 


Yellow 


1 


3 


1 


Red 



Written as a polynomial, f% has the following 
form 

4xjxl + 4xjx|. + x\ + 4x^x|+ 

j k j k 

This polynomial form can then be used to com- 
pute the steady state configuration by solving a 
system of polynomial equations. 
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Figure 3: Dynamics where all cells evolve according to the rules. The center cell determines the 
final outcome. 



6 ODEs with acyclic depen- 
dency graph 

As an aside note we will mention that our main 
result is also valid for dynamical systems with 
differential equations. We now consider ODEs 
that have an acyclic dependency graph, that is, 
no feedback loops in their structure. For any 
ODE with acyclic dependency graph we can see 
easily that its adjacency matrix [/] is a strictly 
triangular matrix, i.e., with zeros in the diagonal 
(otherwise, there would be loops in the graph). 

Theorem 6.1. Consider an ODE with acyclic 
dependency graph, then there is a unique steady 
state and it is asymptotically stable. 

To proof this we need the following lemma. 

Lemma 6.2. Let z : [0, oo) — > M be a continuous 
function such that lim^oo z(t) = zq £ M. If y' = 
z — ky, y(0) = yo with k > 0, then lirrn._s.oo y(t) = 
t& (regardless of the value ofyo). 

Proof. First, notice that from y' = z — ky 
we obtain y(t) = yoe~ kt + f e k ^ s ~^ z(s)ds. 
Furthermore, y(t) — ^ = yoe~ kt + 
Jje fc ( a -*)(*(a) - z Q )ds + f* e^-^zods- f . It is 



2M f 



-fct/a_ 



- + 



jt e Hs-t) Z(jds 



ill 

k 



easy to show that um^oo 

lim^oo yoe~ kt = 0; then, we only need to 
show that lim t _ s , 00 J Q * e k< ^ s ~ t \z(s) — zo)ds = 0. 
First, since lim^oo z{t) = zq, it follows 
that z(s) is bounded (let us say by M) and 
limt-^oo sup s& yt t iA\z(s) — zq\} = 0. On the other 
hand we have 

\J^e k ( s - t \z(s)-z )ds\ < 

fj e k ^\z(s)-z \ds + ft e k ^\z{s) - z \ds 

< 2M fj e k( - s ~ t '>ds+ 
sup se[ t t] {\z(s) - z \} ft e k( - s ~^ds 



sup se[ t_ A {\z{s) - zqW 1 \ kt ' 2 

<2M^ + sup se[ht] {\z(s)-z \}l 
The last expression converges to as t 
and this finishes the proof. 



x 



□ 



Corollary 6.3. Consider x : [0, oo) — > W l , 

f : Bi n — > Bl continuous such that lim^oo x(t) = 
x G R n . If y' = f{x) - ky with k > 0, then 
limt-xx, yit) 



k 



Proof. It is enough to consider z{t) = f(x(t)) 
and apply the lemma above. □ 

When modeling biological systems using 
ODEs, it is common for the functions to have 
natural decay; that is, they are of the form 
x 'i = f( x ) ~ ki x i for i = 1, . . . , n. The depen- 
dency graph of such an ODE is the graph with 
nodes {1, . . . ,n} and an edge from i to j if fj 
depends on X{. 



Proof of theorem 6. 1 Consider an ODE with 
natural decay. Without loss of generality, we 
consider that the adjacency matrix of the de- 
pendency graph is of the form 



... 

* ... 

* * ... 

* * * 



That is, fi = fi(xi, . . . ,xi-i) (/» could depend 
on less variables). For i = 1 we have x[ = fx — 
k\xi where f\ is constant. Then, linn-s.oo x\(t) = 
s\, where s% = j^-. For i = 2 we have x' 2 = 
f2(x±) — k2X2] then, by the corollary, we have 
that lim^oo x 2 (t) = s 2 , where s 2 = 
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By induction, if limt 



Ol, . . . ,Si_l) 

then lim^oo 



and x\ 

Xi(t) = 



oo(xi(t), . 

fi(xi, ■ • 
where s,; 



. ,£i-i(i)) = 
i X{— x) kiXi, 







At the end, we have s = (si, . . . , s n ) G M n such 
that linn._s.oo _ (i) = s (regardless of the value of 
x(0)). It follows that s is the unique steady state 
and that it is asymptotically stable. □ 



7 Discussion 

In this paper we have shown that it is fruit- 
ful to study the relationship between the struc- 
ture and the dynamics of discrete dynamical sys- 
tems by looking at the algebraic properties of 
the mapping V — > Ai from the algebra of PDS 
to the algebra of adjacency matrices. With a 
very straightforward proof we have shown that 
dynamics is very simple in the absence of feed- 
back loops. It is worth observing that this result, 



Corollary 4.1, could also have been obtained as 



a consequence of a more general result that says 
that for the existence of more than one fixed 
point, a positive feedback loop is required in the 
dependency graph, and for periodic orbits to ex- 
ist a negative feedback is necessary (but not suffi- 
cient). This result implies that in order to obtain 
more than one fixed point or periodic orbits, it 
is necessary that the dependency graph of the 
system have feedback loops [S]. However, the 
proof we have given here of this same result is 
very simple and stems from a basic property of 
the mapping V — > Ai rather than complicated 
phase space arguments. It emphasizes our be- 
lief that the proper framework for studying the 
relationship between structure and dynamics of 
PDS is the algebra inherent in this mapping. 
Furthermore, we have shown that one can 



use Corollary 4.1 to draw non-obvious conclu- 
sions about a system of interest. This conclusion 



would be very difficult to arrive at through sim- 
ulations, due to the combinatorial complexity of 
the dynamics on a large grid. 
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